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Abstract 

It is suggested that the colorless systems of interacting soft-gluons in large-rapidity- 
gap events are open dynamical complex systems in which self-organized criticality 
and BTW-clusters play an important role. Theoretical arguments and experimen- 
tal evidences supporting such a statistical approach to deep-inelastic scattering are 
presented. 
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I. Interacting soft gluons in the small-x^ region of DIS 



A number of striking phenomena have been observed in recent deep-inelastic electron-proton 
scattering experiments in the small-x^ region. In particular it is seen, that the contribution 
of the gluons dominates^, and that large-rapidity-gap (LRG) events exist!. The latter shows 
that the virtual photons in such processes may encounter colorless objects originating from 
the proton. 

The existence of LRG events in these scattering processes have attracted much attention, 
and there has been much discussioniHi on problems associated with the origin and/or the 
properties of such colorless objects. Reactions in which "exchange" of such objects dominate 
are known in the literaturelHi as "diffractive scattering processes" . While the concepts and 
methods used by different authors are in general very much different from one another, all 
the authors in describing such processes (experimentalists as well as theorists) seem to agree 
on the followingi (see also Refs. ^|^): (a) Interacting soft gluons play a dominating 

role in understanding the phenomena in the small-x^ region in general, and in describing the 
properties of LRG events in particular, (b) Perturbative QCD should be, and can be, used 
to describe the LRG events associated with high transverse-momentum (p_l) jets which have 
been observed at HERaI and at the TevatronB. Such events are, however, rather rare. For 
the description of the bulk of LRG events, concepts and methods beyond the perturbative 
QCD (for example, Pomeron Modelsi based on Regge Phenomenology) are needed. The 
question, whether or how perturbative QCD plays a role in such non-perturbative approaches 
does not have an unique answer. 

In a previous papeiH, we suggested that the observed dominance of interacting soft 
gluonsS and the existence of LRG event^ in the small- a;^ region are closely related to each 
other, and that the interacting soft gluons may form colored and colorless systems — which 
we called "gluon clusters". Such gluon clusters have finite lifetimes which (in the small-x^ 
region) can be of the same order as the interaction time Tint — the time-interval in which 
the virtual photon 7* "sees" the cluster in the sense that it is absorbed by the charged 
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constituents of the latter. In Ref. P] the hfetime of such a gluon-cluster was estimated by 
using uncertainty principle and kinematical considerations — without any dynamical input. 
In analogy to hadron structure function, a quantity F| which we called "the structure 
function of the gluon cluster Cq" was introduced, and then it was set to be a constant — in 
accordance with the purpose of that paper which is to discuss the kinematical aspects of a 
statistical approach. 

After having seen what phase-space considerations can, and cannot do, we decided to 
go one step further, and study the dynamical aspects of the interacting soft-gluons in these 
scattering processes. In doing so, we realized that the system of interacting soft-gluons 
is extremely complicated. It is not only too complicated (at least for us) to take the de- 
tails of local interactions into account (for example by describing the reaction mechanisms 
in terms of Feynman diagrams), but also too complicated to apply well-known concepts 
and methods in conventional equilibrium statistical mechanics. In fact, having the above- 
mentioned empirical facts about LRG events and the basic properties of gluons prescribed 
by the QCD-Lagrangian in mind, we are readily led to the following picture: 

Such a system is an open dynamical system with many degrees of freedom, and it is 
in general far from equilibrium. This is because, once we accept that the colorless object 
(which the virtual photon encounters) is a system of soft gluons whose interactions are not 
negligible, we are also forced to accept that, in such a system, gluons can be emitted and 
absorbed by the members of the system as well as by gluons and/or quarks and antiquarks 
outside the system (we note in particular that, since the gluons are soft, their density in 
space is high, and the distances between the interacting gluons are in general not short, the 
"running-coupling-constant" can be very large). Furthermore, since in general more than 
one gluons can be emitted or absorbed by the members of the system, the system itself can 
remain to be a color-singlet. This means in particular that, in such a system, neither the 
number of gluons nor the energy of the system can be a conserved quantity. 

Do we see comparable open, dynamical, complex systems in Nature? If yes, what are 
the characteristic features of such systems? 
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II. Characteristic features of open dynamical complex systems 



Open dynamical complex systems are not difficult to find in Nature — at least not in the 
macroscopic world! Such systems have been studied, and in particular the following have 
been observed by Bak, Tang and Wiesenfeld (BTW) some time agoi: Open dynamical 
systems with many degrees of freedom may evolve to self-organized critical states which 
lead to fluctuations extending over all length- and time-scales, and that such fluctuations 
manifest themselves in form of spatial and temporal power-law scaling behaviors showing 
properties associated with fractal structure and flicker noise respectively. 



the following: Dynamical systems with local interacting degrees of freedom can evolve into 
self-organized structures of states which are barely stable. A local perturbation of a critical 
state may "propagate" , in the sense that it spreads to (some) nearest neighbors, and than to 
the next-nearest neighbors, and so on in a "domino effect" over all length scales, the size of 
such an "avalanche" can be as large as the entire system. Such a "domino effect" eventually 
terminates after a total time T, having reached a final amount of dissipative energy and 
having effected a total spatial extension S. The quantity S is called by BTW the "size", 
and the quantity T the "lifetime" of the avalanche — named by BTW a "cluster" (hereafter 
referred to as BTW-cluster). As we shall see in more details later on, it is of considerable 
importance to note that a BTW-cluster cannot, and should not be identified with a cluster in 
the usual sense. It is an avalanche, not a static object with a fixed structure which remains 
unchanged until it decays after a time- interval (known as the lifetime in the usual sense). 

It has been shownlE that the distribution (-D5) of the "size" (which is a measure of the 
dissipative energy, S) and the distribution (-Dt) of the lifetime (T) of BTW-clusters in such 
open dynamical systems obey power-laws: 



BTwi and 



many other authorslSI proposed, and demonstrated by numerical simulations, 




Ds{S) ~ 



(1) 



Dt{T) ~ T 



•—u 



(2) 
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where /i and v are positive real constants. In fact, such spatial and temporal power-law 
scaling behaviors can be, and have been, considered as the universal signals — the "fin- 
gerprints" — of the locally perturbed self-organized critical states in such systems. It is 
expectediS that the general concept of self-organized criticality (SOC), which is comple- 
mentary to chaos, may be the underlying concept for temporal and spatial scaling in a wide 
class of o'pen n on- equilibrium systems — although it is not yet known how the exponents in 
such power law can be calculated analytically. 

SOC has been observed in a large number of open dynamical complex systems in non- 
equilibriumi'0'&^l among which the following examples are of particular interest, because 
they illuminate several aspects of SOC which are relevant for the discussion in this paper. 

First, the well known Gutenberg- Richter law§'§ for earthquakes special case of 
Eq.(l): In this case, S stands for the released energy (the magnitude) of the earthquakes. 
Ds{S) is the number of earthquakes at which an energy S is released. Such a simple law 
is known to be valid for all earthquakes, large (up to 8 or 9 in Richter scale) or small! We 
note, the power-law behavior given by the Gutenberg-Richter law implies in particular the 
following. The question "How large is a typical earthquake?" does not make sense! 

Second, the sandpile experiment JhE which show the simple regularities mentioned in 
Eqs.(l) and (2): In this example, we see how local perturbation can be caused by the addition 
of one grain of sand (note that we are dealing with an open system!). Here, we can also see 
how the propagation of perturbation in form of "domino effect" takes place, and develops 
into avalanches of all possible sizes and durations. The size- and duration-distributions are 
given by Eqs.(l) and (2) respectively. This example is indeed a very attractive one, not only 
because such experiments can be, and have been performed in but also because they 

can be readily simulated on a PcSllsl. 

Furthermore, it has been pointed out, and demonstrated by simple mo dels0'M, that 
the concept of SOC can also be applied to Biological Sciences. It is amazing to see how 
phenomena as complicated as Life and Evolution can be simulated by simple models such 
as the "Game of Life"i and the "Evolution Model" 



Having seen that systems of interacting soft-gluons are open dynamical complex systems, 
and that a wide class of open systems with many degrees of freedom in the macroscopic world 
evolve to self-organized critical states which lead to fluctuations extending over all length- 
and time-scales, it seems natural to ask the following: Can such states and such fluctuations 
also exist in the microscopic world — on the level of quarks and gluons? 



III. Are gluon-clusters hadron-like? 

How can we find out whether the general concept of self-organized criticality (mentioned 
in Section II) plays a role in difFractive deep-inelastic lepton-hadron scattering processes 
(discussed in Section I)? A simple and effective way of doing this, is to check whether 
the "fingerprints" mentioned in Eqs.( |l|) and ( H), which can be considered as the necessary 
conditions for the existence of self-organized criticality, show up in the relevant experiments. 
For such a comparison, we need the spatial and the temporal distributions of the gluon- 
clusters. Hence, an important step in our quantitative study is to obtain these distributions 
directly from the experimental data — if possible, without any theoretical input. Having 
this goal in mind, we now try to express such cluster-distributions in terms of the measured^ 
"diffractive structure function" F^^^\(3,Q'^; xp) = J dtF^^'^\(3,Q'^; xp,t). Here, we note 
that F^^^\(3,Q'^]Xp,t) is relatedH to the differential cross-section for large-rapidity-gap 
events 

■(1 - ir)^2 (/5,<5 ;a;p,^), (3) 



dpdQMxpdt V ^ ■ 2 

in analogy to the relationship between the corresponding quantities [namely d'^a /{dxs dQ'^) 
and F2{xb, Q"^)] for normal deep-inelastic electron-proton scattering events 

'l-y+'-)F2{xB,Q^). (4) 



dxBdQ"^ xbQ^ 2 
The kinematical variables, in particular (3, Q^, xp and xp (in both cases) are directly 
measurable quantities, the definitions of which are shown in Fig.l together with the corre- 
sponding diagrams of the scattering processes. We note that, although these variables are 



Lorentz-invariants, it is sometimes convenient to interpret them in a "fast moving frame", 
for example the electron-proton center-of-mass frame where the proton's 3-momentum P is 
large (i.e. its magnitude \P\ and thus the energy = (|-Pp + M^)^/^ is much larger than 
the proton mass M). While characterizes the virtuality of the space-like photon 7*, xb 
can be interpreted, in such a "fast moving frame" (in the framework of the celebrated parton 
model), as the fraction of proton's energy P° (or longitudinal momentum \P\) carried by 
the struck charged constituent. 

We recall, in the framework of the parton model, F2{x b , Q'^) / x b for "normal events" 
can be interpreted as the sum of the probability densities for the above-mentioned 7* to 
interact with such a charged constituent inside the proton. In analogy to this, the quantity 
F^^^\l3, Q^; xp)/(3 for LRG events can be interpreted as the sum of the probability densities 
for 7* to interact with a charged constituent which carries a fraction (3 = xb/xp oi the energy 
(or longitudinal momentum) of the colorless object, under the condition that the colorless 
object (which we associate with a system of interacting soft gluons) carries a fraction xp 
of proton's energy (or longitudinal momentum). We hereafter denote this charged- neutral 
and color-neutral gluon-system by Cg (in Regge pole modelsi this object is known as the 
"pomeron"). Hence, by comparing Eq. (3) with Eq. (4) and by comparing the two diagrams 
shown in Fig. 1(a) and Fig. 1(b), it is tempting to draw the following conclusions: 

The diffractive process is nothing else but a process in which the virtual photon 7* 
encounters a Cq, and [3 is nothing else but the Bj or ken- variable with respect to Cq (this 
is why it is called xpc in Ref. p[). This means, a diffractive e'p scattering event can be 
envisaged as an event in which the virtual photon 7* collides with "a Cp-target" instead of 
"the proton-target". Furthermore, since Cg is charge-neutral, and a photon can only directly 
interact with an object which has electric charges and/or magnetic moments, it is tempting 
to assign Cg an electromagnetic structure function F2{l3,Q'^), and study the interactions 
between the virtual photon and the quark(s) and antiquark(s) inside Cg. In such a picture 
(which should be formally the same as that of Regge pole models^, if we would replace the 
Cg's by "pomerons") we are confronted with the following two questions: 
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First, is it possible and meaningful to discuss the xp-distributions of the Cq's without 
knowing the intrinsic properties, in particular the electromagnetic structures, of such ob- 
jects? 

Second, are gluon-clusters hadron-like, such that their electromagnetic structures can be 
studied in the same way as those for ordinary hadrons? 

We discuss the second question here, and leave the first question to the next section. We 
note, in order to be able to answer the second question in the affirmative, we need to know 
whether F2^^\l3,Q'^;xp) can be factorized in the form 

F,''^'\l3,Q';xp) = Uxp)F^{[3,Q'). (5) 

Here, fdxp) plays the role of a "kinematical factor" associated with the "target Cq", and 
xp is the fraction of proton's energy (or longitudinal momentum) carried by Cq. [We could 
call fc{xp) "the Cp-flux" — in exactly the same manner as in Regge pole model^, where it 
is called "the pomeron flux".] F2(/5,Q^) is "the electromagnetic structure function of Cq" 
[the counterpart of F2{xb,Q'^) of the proton] which — in analogy to proton (or any other 
hadron) — can be expressed as 

^^^^1^ = E^nQMQ')+mQ% (6) 

where q1{qf) stands for the probability density for 7* to interact with a quark (antiquark) 
of flavor i and electric charge e, which carries a fraction (3 of the energy (or longitudinal 
momentum) of Cq. It is clear that Eq.(6) should be valid for all xp- values in this kinematical 
region, that is, both the right- and the left-hand-side of Eq. (6) should be independent of the 
energy (momentum) carried by the "hadron" Cq. 

Hence, to find out experimentally whether the second question can be answered in the 
affirmative, we only need to check whether the data are in agreement with the assumption 
that F2(/5, Q^) prescribed by Eqs.(5) and (6) exists. For such a test, we take the existing 
datai and plot log[F2^*-^''(/3, Q^; xp)//?] against log/? for different xp-values. We note, under 
the assumption that the factorization shown in Eq.(5) is valid, the /3-dependence for a 



given in such a plot should have exactly the same form as that in the corresponding 
log[F2(/5, vs log/3 plot; and that the latter is the analog of log[F2(xB, (5^)/a;B] vs 

loga;^ plot for normal events. In Fig.2 we show the result of such plots for three fixed Q^- 
values (3.5, 20 and 65 GeV^, as representatives of three different ranges in Q^). Our goal is 
to examine whether or how the /3-dependence of the function given in Eq.(6) changes with 
xp. In principle, if there were enough data points, we should, and we could, do such a plot 
for the data-sets associated with every xp-value. But, unfortunately there are not so much 
data at present. What we can do, however, is to consider the /^-distributions in different 
xp-bins, and to vary the bin-size of xp, so that we can explicitly see whether/how the shapes 
of the /^-distributions change. The results are shown in Fig.2. The /^-distribution in the first 
row, corresponds to the integrated value F-fiP, Q"^) shown in the literatureili. Those in the 
second and in the third row are obtained by considering different bins and/or by varying 
the sizes of the bins. By joining the points associated with a given xp-interval in a plot 
for a given Q^, we obtain the /^-distribution for a Cq carrying approximately the amount 
of energy xpP^, encountered by a photon of virtuality Q^. Taken together with Eq.(6) we 
can then extract the distributions (/?, Q^) and (/5, Q"^) for this Q^-value, provided that 
F2{P,Q'^)/ P is independent of xp. But, as we can see in Fig.2, the existing datJ'i show 
that the xp-dependence of this function is far from being negligible! Note in particular 
that according to Eq.(§), by choosing a suitable /p(xp) we can shift the curves for different 
Xp- values in the vertical direction (in this log-log plot); but we can never change the shapes 
of the f3- distributions which are different for different xp- values! 

In order to see, and to realize, the meaning of the xp-dependence of the distributions 
of the charged constituents of Cg expressed in terms of F2{P,Q'^)/ P in LRG events [see 
Eqs.(5) and (6)], let us, for a moment, consider normal deep-inelastic scattering events 
in the xp-region where quarks dominate (xp > 0.1, say). Here we can plot the data for 
log[F2(xp, (5^)/xp] as a function of logxp obtained at different incident energies (P^'s) 
of the proton. Suppose we see, that at a given Q^, the data for xp-distributions taken 
at different values of P° are very much different. Would it still be possible to introduce 
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F2{xb,Q^) as "the electromagnetic structure function" of the proton, from which we can 
extract the x^-distribution of the quarks qi{xB, Q^) at a given Q^? 

IV. Distributions of the gluon-clusters 

After having seen that the existing data are not in agreement with the picture in which 
the colorless gluon-clusters (cq's) are hadron-like, we now come back to the first question in 
Section III, and try to find out whether it is never-the-less possible and meaningful to talk 
about the xp-distribution of Cq. We shall see in this section, the answer to this question is 
Yes! Furthermore, we shall also see, in order to answer this question in the affirmative, we 
do not need the factorization mentioned in Eq.(5); and we do not need to know whether 
the gluon-clusters are hadron-hke. But, as we shall show later on, it is of considerable 
importance to discuss the second question in understanding the nature of the Cq's. 

In view of the fact that we do use the concept "distributions of gluons" in deep-inelastic 
lepton-hadron scattering, although the gluons do not directly interact with the virtual pho- 
tons, we shall try to introduce the notion "distribution of gluon-clusters" in a similar manner. 
In order to see what we should do for the introduction of such distributions, let us recall the 
following: 

For normal deep-inelastic e'p collision events, the structure function F2{xb,Q^) can be 
expressed in term of the distributions of partons, where the partons are not only quarks 
and antiquarks, but also gluons which can contribute to the structure function by quark- 
antiquark pair creation and annihilation. In fact, in order to satisfy energy-momentum- 
conservation (in the electron- proton system), the contribution of the gluons Xgg{xg, Q^) has 
to be taken into account in the energy-momentum sum rule for all measured Q^-values. Here, 
we denote by g{xg, Q"^) the probability density for the virtual photon 7* (with virtuality Q"^) 
to meet a gluon which carries the energy (momentum) fraction Xg of the proton, analogous 
to qi{xB, Q^) [or qi{xB, Q'^)] which stands for the probability density for this 7* to interact 
with a quark (or an antiquark) of flavor i and electric charge Cj which carries the energy 
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(momentum) fraction xb of the proton. We note, while both xb and Xg stand for energy 
(or longitudinal momentum) fractions carried by partons, the former can be, but the latter 
cannot be directly measured. 

Having these, in particular the energy-momentum sum rule in mind, we immediately see 
the following: In a given kinematical region in which the contributions of only one category 
of partons (for example quarks for xb > 0.1 or gluons for xb < 10^^) dominate, the structure 
function F2{xb, Q^) can approximately be related to the distributions of that particular kind 
of partons in a very simply manner. In fact, the expressions below can be, and have been, 
interpreted as the probability-densities for the virtual photon 7* (with virtuality Q^) to meet 
a quark or a gluon which carries the energy (momentum) fraction Xb ot Xg respectively. 

F2{xb,Q'^) _ 2 / r^2\ F2{Xb,Q'^) , (ry\ 

~ Z^e^ gi(xB,Q ) or ^g{xg,Q). (7) 

X,B j 

The relationship between qi{xB,Q'^), g{xg,Q'^) and F2{xb,Q'^) as they stand in Eq.(^ are 
general and formal (this is the case especially for that between g and F2) in the following 
sense: Both qi{xB, Q"^) and g{xg, Q"^) contribute to the energy-momentum sum rule and both 
of them are in accordance with the assumption that partons of a given category (quarks or 
gluons) dominate a given kinematical region (here xb > 0.1 and xb < 10^^ respectively). 
But, neither the dynamics which leads to the observed Q^-dependence nor the relationship 
between Xg and xb are given. This means, without further theoretical inputs, the simple 
expression for gi^Xg^Q"^) as given by Eq.(7) is practically uselessl 

Having learned this, we now discuss what happens if we assume, in diffractive lepton- 
nucleon scattering, the colorless gluon-clusters (cq's) dominate the small-xs region {xb < 
10^^, say). In this simple picture, we are assuming that the following is approximately true: 
The gluons in this region appear predominately in form of gluon clusters. The interaction 
between the struck Cg and the rest of the proton can be neglected during the 7-Co collision 
such that we can apply impuls-approximation to the Cq's in this kinematical region. That 
is, here we can introduce — in the same manner as we do for other partons (see Eq.^), 
a probability density Dsixpl^jQ"^) for 7* in the diffractive scattering process to "meet" a 
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Cq which carries the fraction xp of the proton's energy P° = (|Pp + M^)^/^ ^ \P\ (where 
P is the momentum and M is the mass of the proton). In other words, in diffractive 
scattering events for processes in the kinematical region xp < 10^^, we should have, instead 
of g{xg,Q^), the following: 

F,''^'\p,Q^;xp) 



Xp 



Ds{xp\(3,Q'). (8) 



Here, xpP^ is the energy carried by Cq, and j3 indicates the corresponding fraction carried by 
the struck charged constituent in Cq. In connection with the similarities and the differences 
between qi{xp, Q^), g^xp, Q^) in (0) and Ds{xp\P, Q"^) in (H), it is useful to note in particular 
the significant difference between Xg and xp, and thus that between the x^-distribution 
g{xg, Q"^) of the gluons and the xp-distribution Ds{xp\P, Q"^) of the Cq's: Both Xg and xp are 
energy (longitudinal momentum) fractions of charge-neutral objects, with which 7* cannot 
directly interact. But, in contrast to x^, xp can be directly measured in experiments, namely 
by making use of the kinematical relation 

and by measuring the quantities Q^, and in every collision event. Here, Q, and 
W stand respectively for the invariant momentum-transfer from the incident electron, the 
invariant-mass of the final hadronic state after the 7* — Cq collision, and the invariant mass 
of the entire hadronic system in the collision between 7* and the proton. Note that xp = 
(3xp, hence jS is also measurable. This means, in sharp contrast to g{xg,Q'^), experimental 
information on Ds{xp\[3,Q'^) in particular its xp-dependence can be obtained — without 
further theoretical inputsl 



V. The first SOC-fingerprint: Spatial scaling 

We mentioned at the beginning of Section HI, that in order to find out whether the concept 
of SOC indeed plays a role in diffractive DIS we need to check the fingerprints of SOC 
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shown in Section II, and that such tests can be made by examing the corresponding cluster- 
distributions obtained from experimental data. We are now ready to do this, because we have 
learned in Sections III and IV, that it is not only meaningful but also possible to extract xp- 
distributions from the measured diffractive structure functions, although the gluon-clusters 
cannot be treated as hadrons. In fact, as we can explicitly see in Eqs.(8) and (9), in order 
to extract the xp-dependence of the gluon-clusters from the data, detailed knowledge about 
the intrinsic structure of the clusters are not necessary. 

Having these in mind, we now consider Ds as a function of xp for given values of j3 
and Q^, and plot F^^^\j3, Q^; xp)/xp against xp for different sets of (3 and Q^. The results 
of such log-log plots are shown in Fig. 3. As we can see, the data! suggest that the 
probability-density for the virtual photon 7* to meet a color-neutral and charged-neutral 
object Cq with energy (longitudinal momentum) fraction xp has a power-law behavior in 
Xp, and the exponent of this power- law depends very little on and (3. This is to be 
compared with Ds{S) in Eq.( [1|), where S, the dissipative energy (the size of the BTW- 
cluster) corresponds to the energy of the system Cq. The latter is xpP^, where P° is the 
total energy of the proton. 

It means, the existing datai show that Ds{xp\P,Q'^) exhibits the same kind of power- 
law behavior as the size-distribution of BTW-clusters. This result is in accordance with 
the expectation that self-organized critical phenomena may exist in the colorless systems of 
interacting soft gluons in diffractive deep-inelastic electron-proton scattering processes. 

We note, up to now, we have only argued (in Section I) that such gluon-systems are 
open, dynamical, complex systems in which SOC may occur, and we have mentioned (in 
Section II) the ubiquitousness of SOC in Nature. Having seen the first piece of experimental 
evidence that one of the necessary conditions for the existence of SOC is satisfied, let us now 
take a second look at the colorless gluon-systems from a theoretical point of view: Viewed 
from a "fast moving frame" which can for example be the electron-proton c.m.s. frame, such 
colorless systems of interacting soft gluons are part of the proton (although, as color-singlets, 
they can also be outside the confinement region). Soft gluons can be intermittently emitted 
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or absorbed by gluons in such a system, as well as by gluons, quarks and antiquarks outside 
the system. The emission- and absorption-processes are due to local interactions prescribed 
by the well-known QCD-Lagrangian (here "the running coupling constants" are in general 
large, because the distances between the interacting colored objects cannot be considered 
as "short"; remember that the spatial dimension of be much larger than that of 

a hadron!). In this connection, it is however very useful to keep in mind that, due to the 
complexity of the system, details about the local interactions may be relatively unimportant, 
while general and/or global features — for example energy-flow between different parts 
(neighbors and neighbor's neighbors . . .) of the system — may play an important role. 
How far can one go in neglecting dynamical details when one deals with such open 



in modelling some of the essential aspects of The Evolution in Nature. They consider 
the "fitness" of different "species", related to one another through a "food chain", and 
assumed that the species with the lowest fitness is most likely to disappear or mutate at 
the next time-step in their computer simulations. The crucial step in their simulations 
that drives evolution is the adaption of the individual species to its present environment 
(neighborhood) through mutation and selection of a fitter variant. Other interacting species 
form part of the environment. This means, the neighbors will be influenced by every time- 
step. The result these authors obtained strongly suggests that the process of evolution is a 
self-organized critical phenomenon. One of the essential simplifications they made in their 
evolution mo del^l^'i is the following: Instead of the explicit connection between the fitness 
and the configuration of the genetic codes, they use random numbers for the fitness of the 
species. Furthermore, as they have pointed out in their papers, they could in principle have 
chosen to model evolution on a less coarse-grained scale by considering mutations at the 
individual level rather than on the level of species, but that would make the computation 
prohibitively difficult. 

Having these in mind, we are naturally led to the questions: Can we consider the creation 
and annihilation processes of colorless systems of interacting soft gluons associated with a 



complex systems? In order to see this, let us recall 




succeeded 
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proton as "evolution" in a microscopic world? Before we try to build models for a quan- 
titative description of the data, can we simply apply the existing evolution model^llll to 
such open, dynamical, complex systems of interacting soft-gluons, and check whether some 
of the essential features of such systems can be reproduced? 

To answer these questions, we now report on the result of our first trial in this direc- 
tion: Based on the fact that we know very little about the detailed reaction mechanisms in 
such gluon-systems and practically nothing about their structures, we simply ignor them, 
and assume that they are self-similar in space (this means, colorless gluon-clusters can be 
considered as clusters of colorless gluon-clusters and so on). Next, we divide them in an 
arbitrary given number of subsystems Sj (which may or may not have the same size). Such 
a system is open, in the sense that neither its energy e^, nor its gluon- number has a 
fixed value. Since we do not know, in particular, how large the e^'s are, we use random 
numbers. As far the nj's are concerned, since we do not know how these numbers are asso- 
ciated with the energies in the subsystems Sj, except that they are not conserved quantities, 
we just ignor them, and consider only the e^'s. As in Ref . ||14]| or in Ref. |T5|, the random 



number of this subsystem as well as those of the fixed0 or random (see the first paper of 
Ref. ||15||) neighbors will be changed at every time-step. Note, this is how we simulate the 
processes of energy flow due to exchange of gluons between the subsystems, as well as those 
with gluons/quarks/antiquarks outside the system. In other words, in the spirit of Bak 
and Sneppen0 we neglecting the dynamical details totally. Having in mind that, in such 
systems, the gluons as well as the subsystems (sj's) of gluons are virtual (space-like), we 
can ask: "How long can such a colorless subsystem Sj of interacting soft gluons exist, which 
carries energy e^?" According to the uncertainty principle, the answer should be: "The 
time interval in which the subsystem Si can exist is proportional to l/si, and this quantity 
can be considered as the lifetime Tj of Sj." In this sense, the subsystems of colorless gluons 
are expected to have larger probabilities to mutate because they are associated with higher 
energies, and thus shorter "lifetimes". Note that the basic local interaction in this self- 
organized evolution process is the emission (or absorption) of gluons by gluons prescribed 
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by the QCD-Lagrangian — although the detailed mechanisms (which can in principle be 
explicitly written down by using the QCD-Lagrangian) do not play a significant role. 

In terms of the evolution model0lli we now call Si the "species" and identify the corre- 
sponding lifetime Tj as the "fitness of s" . Because of the one-to-one correspondence between 
Ti and Ei, where the latter is a random number, we can also directly assign random numbers 
to the Tj's instead. From now we can adopt the evolution mo del0@ and note that, at the 
start of such a process (a simulation), the fitness on average grow, because the least fit are 
always eliminated. Eventually the fitness do not grow any further on average. All gluons 
have a fitness above some threshold. At the next step, the least fit species (i.e. the most 
energetic subsystem Sj of interacting soft gluons), which would be right at the threshold, 
will be "replaced" and starts an avalanche (or punctuation of mutation events), which is 
causally connected with this triggering "replacement". After a while, the avalanche will 
stop, when all the fitnesses again will be over that threshold. In this sense, the evolution 



goes on, and on, and on. As in Refs.|]n[ and [|T^], we can monitor the duration of every 
avalanche, that is the total number of mutation events in everyone of them, and count how 
many avalanches of each size are observed. The avalanches mentioned here are special cases 
of those discussed in Section II. Their size- and lifetime-distributions are given by Eq.(l) 
and Eq.(2) respectively. Note in particular that the avalanches in the Bak-Sneppen model 
correspond to sets of subsystems Sj, the energies (ej) of which are too high "to be fit for the 
colorless systems of low-energy gluons". It means, in the proposed picture, what the virtual 
photon in deep-inelastic electron-proton scattering "meet" are those "less fit" one — those 
who carry "too much" energy. In a geometrical picture this means, it is more probable for 
such "relatively energetic" colorless gluons-clusters to be spatially further away from the 
(confinement region of) the proton. 

There exists, in the mean time, already several versions of evolution mo dels0§ based on 
the original idea of Bak and Sneppen0 Although SOC phenomena have been observed in all 
these cases0'0'@, the slopes of the power-law distributions for the avalanches are different 
in different models — depending on the rules applied to the mutations. The values range 
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from approximately —1 to approximately —2. Furthermore, these modelstJallj seem to 
show that neither the size nor the dimension of the system used for the computer simulation 
plays a significant role. 

Hence, if we identify the colorless charge-neutral object Cq encountered by the virtual 
photon 7* with such an avalanche, we are identifying the lifetime of Cq with T, and the "size" 
(that is the total amount of dissipative energy in this "avalanche") with the total amount 
of energy of Cq. Note that the latter is nothing else but xpP^, where P° is the total energy 
of the proton. This is how and why the ^-distribution in Eq. ([^) and the xp-distribution of 
Ds{xp\i3,Q^) in Eq.® are related to each other. 

VI. The second fingerprint: Temporal scaling 

In this section we discuss in more detail the effects associated with the time-degree-of- 
freedom. In connection with the two questions raised in Section III, one may wish to know 
why the parton-picture does not always work when we apply it in a straightforward manner 
— not only to hadrons but also to gluon-clusters. The answer is very simple: The time-degree 
of freedom cannot be ignored when we wish to find out whether impulse-approximation is 
applicable, and the applicability of the latter is the basis of the parton-model. We recall 
that, when we apply this model to stable hadrons, the quarks, antiquarks and gluons are 
considered as free and stable objects, while the virtual photon 7* is associated with a given 
interaction-time Tint(Q^,a;s) characterized by the values and xb of such scattering pro- 
cesses. We note however that, this is possible only when the interaction-time Tint is much 
shorter than the corresponding time-scales (in particular the average propagation-time of 
color-interactions in hadron). Having these in mind, we see that, we are confronted with the 
following questions when we deal with gluon-clusters associated with finite lifetimes: Can 
we consider the Cq's as "/ree" and "sto6/e" particles when their lifetimes are shorter than 
the interaction-time Tint(<5^, xp)? Can we say that a 7* — Cg collision process takes place, 
in which the incident 7* is absorbed by one a or a system of the charged constituents of Cq, 
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when the hfetime T of Cq is shorter ihasi Tint(Q^, a^s)? 

Since the notion "stable objects" or "unstable objects" depends on the scale which is 
used in the measurement, the question whether a Cq can be considered as a parton (in the 
sense that it can be considered as a free "stable object" during the 7*-Cq interaction) depends 
very much on on the interaction-time TintiQ"^, xb)- Here, for given values of Q^, xb-, and thus 
TintiQ^ , X b) , only those Cq's whose lifetime (T's) are greater than TintiQ"^ iXb) can absorb 
the corresponding 7*. That is to say, when we consider diffractive electron- proton scattering 
in kinematical regions in which Cq's dominate, we must keep in mind that the measured 
cross-sections (and thus the diffractive structure function F^^^^) only include contributions 
from collision-events in which the condition T > Tint{Q^,XB) is satisfied! 

We note that Tint can be estimated by making use of the uncertainty principle. In fact, 
by calculating in the above-mentioned reference frame, we obtain 

which implies that, for given \P\ and values. 

Tint oc Xb, for Xb < 1- (11) 

This means, for diffractive e~p scattering events in the small-a;^ region at given \P\ and 
values, xb is directly proportional to the interaction time Tint- Taken together with the 
relationship between Tint and the minimum lifetime T(mim) of the Cq's mentioned above, 
we reach the following conclusion: The distribution of this minimum value, T(min) of the 
Cq's which dominate the small-x^ {xb < 10~^, say) region can be obtained by examining 
the xe-dependence of F^^^\/3,Q'^]Xp)//3 discussed in Eqs. (5), (6) and in Fig. 2. This is 
because, due to the fact that this function is proportional to the quark (antiquark) distribu- 
tions qf{qi'^) which can be directly probed by the incident virtual photon 7*, by measuring 
F^^^\f3, Q^,xp)//3 as a, function of xb = /3xp, we are in fact asking the following questions: 
Do the distributions of the charged constituents of Cq depend on the interaction time Tint, 
and thus on the minimum lifetime T(min) of the to be detected gluon-clusters ? We use 
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the identity xb = f3xp and plot the quantity F^^^\(3, Q^; xp)/ (3 against the variable xb for 
fixed values of /3 and Q^. The result of such a log- log plot is given in Fig.4. It shows not 
only how the dependence on the time-degree-of-freedom can be extracted from the existing 
datai, but also that, for all the measured values of (3 and Q^, the quantity 

is approximately independent of /5, and independent of . This strongly suggests that 
the quantity given in Eq.(|12|) is associated with some global features of Cq — consistent 
with the observation made in Section III which shows that it cannot be used to describe the 
structure of Cq. This piece of empirical fact can be expressed by setting p{xb\I3i Q"^) ~ p{xb)- 
By taking a closer look at this log-log plot, as well as the corresponding plots for different 
sets of fixed (3- and Q^-values (such plots are not shown here, they are similar to those in 
Fig. 3), we see that they are straight lines indicating that p{xb) obeys a power-law. What 
does this piece of experimental fact tell us? What can we learn from the distribution of the 
lower limit of the lifetimes (of the gluon-systems Cq's)? 

In order to answer these questions, let us, for a moment, assume that we know the 
lifetime-distribution Dt(T) of the Cq's. In such a case, we can readily evaluate the integral 



I[n^,{xB)] = / DT{T)dT, (13) 

and thus obtain the number density of all those clusters which live longer than the interaction 
time Ti^ti^B)- Hence, under the statistical assumption that the chance for a 7* to be absorbed 
by one of those Cg's of lifetime T is proportional to Dt{T) (provided that T\iit{Q'^ ■, x b) < 
otherwise this chance is zero), we can then interpret the integral in Eq.(13) as follows: 
-^[Tint(Q^, Xb)] oc p{xb) is the probability density for 7* [associated with the interaction-time 
Tmt(a^B)] to be absorbed by Cq's. Hence, 

Drixs) oc -^p{xb)- (14) 

dXB 

This means in particular, the fact that p{xb) obeys a power-law in implies that Dt{T) 
obeys a power-law in T. Such a behavior is similar to that shown in Eq.( 0). In order to see 



the quality of this power-law behavior of Dt, and the quality of its independence of and 
/3, we compare the above-mentioned behavior with the existing datal. In Fig. 5, we show 
the log-log plots of d/dxB[p{xB)] against xb- We note that d/dxB[p{xB)] is approximately 
F^^^\(3, xbI I3)/{I3xb)- The quality of the power-law behavior of Dt is exphcitly shown 
in Fig. 5. 

VII. Q^-dependent exponents in the power-laws? 

We have seen, in Sections V and VI, that in diffractive deep-inelastic electron-proton scat- 
tering, the size- and the lifetime-distributions of the gluon- clusters obey power-laws, and 
that the exponents depend very little on the variables (3 and . We interpreted the power- 
law behaviors as the fingerprints of SOC in the formation processes of such clusters. Can 
such approximately independence (or weak dependence) of the exponents on Q"^ and (3 be 
understood in a physical picture based on SOC? In particular, what do we expect to see in 
photoproduction processes where the associated value for is zero? 

In order to answer these questions, let us recall the space-time aspects of the collision 
processes which are closely related to the above-mentioned power-law behaviors. Viewed 
in a fast moving frame (e.g. the c.m.s. of the colliding electron and proton), the states 
of the interacting soft gluons originating from the proton are self-organized. The colorless 
gluon-clusters caused by local perturbations and developed through "domino effects" are 
BTW-clusters. That is, they are avalanches (see Sections I and V), the size-distribution of 
which [see Eqs.(8) and (1)] are given by Fig.3. This explicitly shows that there are gluon- 
clusters of all sizes, because a power-law size-distribution implies that there is no scale in size. 
Recall that, since such clusters are color-singlets, their spatial extensions can be much larger 
than that of the proton, and thus they can be "seen" also outside the proton by a virtual 
photon originating from the electron. In other words, what the virtual photon encounters is 
a cloud of colorless gluon-clusters spatially extended in- and outside the proton. 

The virtual photon, when it encounters a colorless gluon-cluster, will be absorbed by 
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the charged constituents (quarks and antiquarks due to fluctuation of the gluons) of the 
gluon-system. Here it is useful to recall that in such a space-time picture, is inversely 
proportional to the transverse size, and is a measure of the interaction time [See Eqs. 
(10) and (11) in Section VI] of the virtual photon. It is conceivable, that the values for the 
cross-sections for virtual photons (associated with a given and a given xb) to collide with 
gluon-clusters (of a given size and a given lifetime) may depend on these variables. But, since 
the processes of self-organization (which produce such gluon-clusters) take place independent 
of the virtual photon (which originates from the incident electron and enters "the cloud" to 
look for suitable partners), the power- law behaviors of the size- and lifetime- distributions 
of the gluon-clusters are expected to be independent of the properties associated with the 
virtual photon. This means, by using 7*'s associated with different values of to detect 
clusters of various sizes, we are moving up or down on the straight lines in the log-log plots 
for the size- and lifetime distributions, the slopes of which do not change. In other words, the 
approximative Q^-independence of the slope is a natural consequence of the SOC picture. 

As far as the /3-dependence is concerned, we recall the results obtained in Sections III 
and IV, which exphcitly show the following: The gluon-clusters (cq's) can not be considered 
as hadrons. In particular, it is neither possible nor meaningful to talk about "the electro- 
magnetic structure of the gluon-cluster" . This suggests, by studying the /3-dependence of 
the "diffractive structure functions" we cannot expect to gain further information about the 
structure of the gluon-clusters or further insight about the reaction mechanisms. 

Having seen these, we try to look for measurable quantities in which the integrations over 
/? have already been carried out. A suitable candidate for this purpose is the differential 
cross-section 



Together with Eqs. (3) and (8), we see that this cross-section is nothing else but the effective 
/9-weighted xp-distribution Ds{xp\C^ ^(S) of the gluon-clusters. Note that the weighting 




(15) 
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factors shown on the right-hand-side of Eq.(15) are simply results of QED! Next, we use the 
datai for F^''^^ which are available at present, to do a log-log plot for the integrand of the 
expression in Eq.(15) as a function of xp for different values of (3 and Q^. This is shown 
in in Fig. 6a. Since the absolute values of this quantity depend very much, but the slope of 
the curves very little on f], we carry out the integration as follows: We first fit every set 
of the data separately. Having obtained the slopes and the intersection points, we use the 
obtained fits to perform the integration over f3. The results are shown in the 
/ 1 \ 

plots of Fig. 6b. These results show the Q^-dependence of the slopes is practically negligible, 
and that the slope is approximately —1.95 for all values of Q^. 

Furthermore, in order to see whether the quantity introduced in Eq.(15) is indeed useful, 
and in order to perform a decisive test of the Q^-independence of the slope in the power-law 
behavior of the above-mentioned size-distributions, we now compare the results in deep- 
inelastic scattering with those obtained in photoproduction0, where LRG events have also 
be observed. This means, as in diffractive deep-inelastic scattering, we again associate the 
observed effects with colorless objects which are interpreted as system of interacting soft 
gluons originating from the proton. In order to find out whether it is the same kind of 
gluon-clusters as in deep-inelastic scattering, and whether they "look" very much different 
when we probe them with real {Q^ = 0) photons, we replot the existing da/dM^ data@ for 
photopro duct ion experiments performed at different total energies, and note the kinematical 
relationship between M|-, W'^ and xp for Q"^ <C and |t| <C Mj.: 
Ml+t Ml 

The result of the corresponding 

log ( r~~l versus log(a;p) 

\xp dxp J 

plot is shown in Fig. 7. The slope obtained from a least-square fit to the existing data0 is 
-1.98 ±0.07. 
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The results obtained in diffractive deep-inelastic electron-proton scattering and that for 
diffractive photoproduction strongly suggest the following: The formation processes of gluon- 
clusters in the proton is due to self-organized criticality, and thus the spatial distributions 
of such clusters — represented by the xp-distribution — obey power-laws. The exponents 
of such power-laws are independent of Q^. Since can be interpreted, in a geometrical 

picture, as a measure for the transverse size of the incident virtual photon, the observed 
Q^-independence of the exponents can be considered as further evidence for SOC — in the 
sense that the self-organized gluon-cluster formation processes take place independent of the 
virtual photon (which is "sent in" to detect the clusters). 

VIII Concluding remarks 

The existence of large rapidity gap (LRG) event^H in deep-inelastic electron-proton scatter- 
ing is one of the most striking features, if not the most striking feature of the experimental 
data obtained in the small-xs [xb < 10~^, say) region. Taken together with the empirical 
that gluons dominate in this kinematical region and that their interactions are not 
negligible, it seems quite natural to think, that such events are due to collisions between the 
virtual photons originated from the lepton and colorless gluon-systems originating from the 
proton. 

What we propose in the present paper is a statistical approach to study such colorless 
gluon-systems. The reasons, why we think such an approach is useful, can be summarized 
as follows: 

First, a number of theoretical arguments and experimental indications suggest that such 
a system of interacting soft-gluons is a system with the following properties: (a) It is a 
complex system with many degrees of freedom, because in general it has a large — unknown 
— number of gluons. (b) It is an open system. This is because the members of a colorless 
gluon-system may interact (through emission and/or absorption of soft gluons) not only 
with one another, but also with gluons and/or quarks and antiquarks outside the system. 
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Thus, due to such interactions, neither the gluon-number nor the energy of this system can 
remain constant, (c) It is neither in chemical nor in thermal equilibrium. This is because, 
as we can for example see in the analysis shown in Section III, it is not possible to consider 
the colorless gluon-cluster Cg as a hadron-like object which has a given structure. In this 
sense, we are forced to consider it as a dynamical system — probably very far from thermal 
and chemical equilibria, (d) The basic interactions between the members of the system, as 
well as those between a member and quarks or gluons outside the system, are local. In fact, 
they are explicitly given by the well-known QCD-Lagrangian. But, as it is often the case in 
complex systems, whether the local dynamical details or the general global features of the 
system plays a more significant role is a different question. 

Second, it has been proposed by Bak, Tang and Wiesenfeldi§ some time ago, that a 
wide class of open dynamical complex systems far from equilibrium may evolve in a self- 
organized manner to critical states, which give rise to spatial and temporal power-law scaling 
behaviors. Such scaling behaviors are universal and robust, in fact they can be considered 
as the "fingerprints" of self-organized criticality (SOC). In the macroscopic world, there 
are many open dynamical complex systems which show this kind of scaling behavior^0. 
Under the condition (see above) that the colorless system of interacting gluons can indeed 
be considered as an open, dynamical, complex system, it would be of considerable interest 
to see whether there can be self-organized criticality also in the microscopic world — at the 
level of gluons and quarks. 

Third, by using the existing data for deep inelastic electron-proton scattering! and those 
for photoproductionEi, where colorless systems of interacting soft-gluons are expected to 
play a dominating role, we checked the above-mentioned fingerprints. The obtained results 
show that the above-mentioned characteristic features for SOC indeed exist. Furthermore, 
it is seen that the relevant exponents in such power-laws are the same for different reac- 
tions. The existence of SOC in systems of interacting soft gluons in such reactions has a 
number of consequences. It seems worthwhile to study them in more detail. In particular, it 
would be very helpful to build realistic models and/or cellular automata to do quantitative 
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calculations. 

Fourth, based on the obtained results in particular the validity of the power-law behav- 
iors, the physical picture for a colorless gluon-cluster should be as follows: It is not a hadron 
with a given structure. It has neither a typical size, nor a typical lifetime, and its structure 
is changing all the time. In fact, it has much in common with an earthquake or an avalanche 
(mentioned in more detail in Sections II, IV and V). Can we learn more about these ob- 
jects by studying other reactions? Can we use the same concepts and methods to treat 
hadron-hadron and hadron- nucleus collision processes? It is known that "the exchange of 
colorless objects" plays an important role also in diffractive hadron-hadron collisions. Shall 
we see this kind of power-law behaviors also in diffractive inelastic hadron-hadron scattering 
processes? Studies along this line are in progress. The results will be pubhshed elsewhere, 
when they are ready. 
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FIGURES 

Fig. 1. The well-known Feynman diagrams (a) for diffractive and (b) for normal deep-inelastic 
electron-proton scattering are shown together with the relevant kinematical variables which describe 
such processes. 

Fig. 2. F2^^\(3,Q^;xp)/(3 is plotted as a function of /? for given xp-intervals and for fixed 
Q^-values. The data are taken from Ref.[^. The lines are only to guide the eye. 

Fig. 3. F^^^\p, Q^; xp)/xp is plotted as a function of xp for different values of fi and Q^. The 
data are taken from Ref.|p. 

Fig. 4. F^^^\i3,Q'^; Xp)/ (3 is plotted as a function of in the indicated /3- and Q^-ranges. The 
data are taken from Ref.j^. 

Fig. 5. F^^^\f3,Q'^;xB/P)/{(3xB) is plotted as a function of xp for fixed /3- and Q^-values. The 
data are taken from Ref.|^. 

Fig. 6. Figure 6: (a) {l/xp)d^a^/dpdQ'^dxp is plotted as a function of xp in different bins of (3 
and Q^. The data are taken from Ref.[p. (b) {l/xp)d'^a^ /dQ'^dxp is plotted as a function of xp 
in different bins of . The data are taken from Ref. . 

Fig. 7. Figure 7: {l/xp)da/dxp for photoproduction j + p^X + p is plotted as a function of 
Xp. The data are taken from Ref.||l^. Note that the data in the second paper are given in terms of 
relative cross sections. Note also that the slopes of the straight-lines are the same. The two dashed 
lines indicate the lower and the upper limits of the results obtained by multiplying the lower solid 
line by ct^qi = 154 it 16(stat.) it 32(syst.);u6. This value is taken from the third paper in Ref.|16|. 
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Fig.1 (b) 
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